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Abstract 

We have employed the Density Functional Theory formalism to investigate the nematic-isotropic 
capillary transitions of a nematogen confined by walls that favor antagonist orientations to the liq- 
uid crystal molecules (hybrid cell) . We analyse the behavior of the capillary transition as a function 
of the fluid-substrate interactions and the pore width. In addition to the usual capillary transi- 
tion between isotropic-like to nematic-like states, we find that this transition can be suppressed 
when one substrate is wet by the isotropic phase and the other by the nematic phase. Under 
this condition the system presents interface-like states which allow to continuously transform the 
nematic-like phase to the isotropic-like phase without undergoing a phase transition. Two different 
mechanisms for the disappearance of the capillary transition are identified. When the director of 
the nematic-like state is homogeneously planar-anchored with respect to the substrates, the cap- 
illary transition ends up in a critical point. This scenario is analogous to the observed in Ising 
models when confined in slit pores with opposing surface fields which have critical wetting transi- 
tions. When the nematic-like state has a linearly distorted director field, the capillary transition 
continuously transforms in a transition between two nematic-like states. 
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I. INTRODUCTION 



The effect on confinement on simple fluids has been vastly investigated in the past. Most 
of studies are related to the effect of the confinement on the vapor-liquid transition, the so- 
called capillary transition. In the case of symmetric walls, where both walls attract or repel 
the molecules with the same strength, the capillary transition in large pores is governed by 
the well-known Kelvin law^. In addition, it may appear surface phenomenology that can 
affect the global phase diagram. Often this surface phenomenology is related to wetting 
states found in the single wall cases. The latter phenomenon is related to the spreading of 
drops of vapor (liquid) on the surface in coexistence with the bulk liquid (vapor). The pres- 
ence of these surface states can be considered in the description of the phase behavior in the 
confined system. The generalized Kelvin law for small pores (where surface phenomenology 
is more dominant) takes into account these effects so it is still possible to determine the 
vapor-liquid transition for small pores^. However, the predictions of the generalized Kelvin 
law reveal that the wetting phenomena does not affect essentially the phase diagram of the 
confined problem so they are often discarded. An interesting system is the hybrid cell (also 
known as asymmetric cell) where one substrate or wall is repulsive and the other attractive 
(competing walls). In this case, the wetting properties present in the single wall cases play a 
crucial role in the determination of the critical point in the confined problem. Here, given a 
pore width, the coexistence of two phases can only occur for temperatures below T\y, which 
corresponds to the total wetting temperature of the single wall casea^"^. 

These studies have elucidated which mechanisms are involved in the surface transitions 
and how the properties of the substrate can determine the phase behavior within the pore. 
With respect to complex fluids like liquid crystal (LC) there is an increasing interest in 
the understanding of the surface phenomena due to its crucial role in the advance of the 
LC technology. For the case of a single wall (semi-infinite problem) and symmetric pores 
favoring an specific orientation on the molecules, theoretical approaches^'^'^'^'Siilliil as well 
as computer simulation a^^i^^i^^i^^i^^ reveal that the interaction wall-fluid particle can affect 
strongly the ordering not only close to the surface but also far from it. Additionally to finite 
size effects, novel phenomena often appear due to the combined presence of anchoring and 
orientational wetting transitions. The interplay between both surface phenomena and its 
dependence on wall properties have been recently studied by us employing a microscopic 
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Density Functional Theory (DFT) approach for a simple model of liquid crystali^. The 
effect of confinement by symmetric walls on these surface transitions was also studiecii^. 
Our findings revealed that as in the case of simple fluids, the finite-size effects for not very 
narrow pores on the phase diagram are not of qualitative significant importance and only 
the location of the capillary transitions are affected by the confinement. However, recent 
theoreticali^^ and simulation results^^ for hard-spherocylinder models show that for very 
narrow symmetrical pores the NI capillary transition may end up in a critical point. 

In a recent paper—, we extended the study to the case with an hybrid cell where the 
surfaces favor antagonist orientations. This geometry has in LCs a particular interest due 
to the possible technological applications. This fact has stimulated an increasing inter- 
est by experimentalists and theoretician o^^i^^i^^ . In particular, different nematic-like states 
with homogeneous and linearly distorted director fields were found to coexist under con- 
finement. In the present paper, our goal is to complete this study with the investigation 
of the nematic-isotropic (NI) transition within the hybrid cell. Quintana et al investigated 
a similar problem employing a Landau-type theory^i. They claimed that the NI transi- 
tion either appears unaffected and its temperature T^i remains unchanged or the transition 
disappears. Basically what determined the existence of the transition was the value of fluid- 
wall interaction strength with respect to the value for the first-order wetting transition for 
the semi-infinite systems. For strength values below this threshold the transition occurs for 
all wall separations. In contrast, above the threshold, the transition disappears for large 
pore widths, where interfacial-like states appear. These states are related to the wetting 
states by isotropic and nematic phase that each wall present in the semi-infinite geometry, 
respectively. 

In this paper we address the same problem of the behavior of the NI capillary transition 
and its possible disappearance from a more microscopic point of view by using a DFT 
approach. Our paper is organized as follows. First we briefly present our liquid crystal model 
in Section II. In Section III we revisit and complete the study of a substrate which favors 
either planar or homeotropic anchoring in contact with a nematogen fluid at its isotropic or 
nematic phase. In particular, the different wetting regimes are found. This phenomenology 
will be a crucial guide to understand the behavior of the fluid under confinement. Section 
IV is devoted to the study of the fluid adsorbed in a hybrid cell, paying special attention to 
the disappearance of the NI capillary transitions and the mechanisms involved in it. Finally 
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we end up with our conclusions. 



II. THE MODEL 

The theoretical model is a standard generalized van der Waals theory based on a pertur- 
bative expansion, using a hard-sphere (HS) fluid as reference system"^. Details on the physical 
basis of the model and how to obtain its solutions numerically can be found elsewher o^^i^^ . 
Our starting point is the grand potential functional per unit system area A, Q[p]/A, whose 
functional minimum with respect to the one-particle distribution function, p(r,0), which 
depends on both molecular positions r and orientations Cl, gives the equilibrium structure of 
the interface. This function, p(r, ft) = p{z)f{z, O), contains a mass distribution p{z) and an 
angular distribution /(2,0). These quantities vary locally with the distance from 2; = to 
z = H, where H is the pore width. The expression for fi[p], in a mean field approximation 
is, 

n[p] = Fr[p] + - [ [ [ [ dvdr'dndn'p{r,n)p{v\n') 
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xvir~v',n,n')~ j J drdnpiv,n)[ix-Vwir,Q)], (1) 

where fi is the chemical potential and 

Fr[p] = J drUM^)) + ksT J drp{r){H4nfiz,n))) (2) 

is the reference system free energy. In the above expression fhs{p{r)) is the hard-sphere free 
energy density of a uniform fluid with a density equal to the local density at r and (...) is an 
angular average. The attractive potential v contains anisotropic (dispersion) forces driving 
the liquid-crystalline behavior of the model material: 

v{r, h, h') = VA{r) + VBir)P2in ■ Q') (3) 



+vc{r) 



P2(0-r) + P2(fi'-r) 



where r = r/r, and r is the center-of mass distance. The functions VA{r), vb{t)-, vc{t) are the 
radial contribution to the potential as a result of considering dispersion-type forces. In this 
work we choose them to have a simple Yukawa form, i.e. fj(r) = — e^exp (— Aj(r — d))/r 
for r > d, and Vi{r) = otherwise, where d is the diameter of a hard sphere. We note 
that a positive value of ec favors the molecule to orientate parallel to the NI interface 



and a negative value favors the perpendicular orientation to NI interface. This anisotropy 
in the orientational configuration of the nematic molecules is essential to understand the 
phenomenology in the problem with the walls. 

The walls are modelled via the following potentials: 

{+00 z <0 

(4) 

where 6 is the angle between the molecule axis and the z-axis normal to the substrate. 
For ejy > 0, this potential favors the molecules to align perpendicularly to the surface 
(homeotropic anchoring). The potential corresponding to the other substrate is: 

{+00 z > H 

(5) 

which favors the molecules to align in a plane parallel to the surface for > 0. However, 
this preferred alignment corresponds to a random planar anchoring, since there is not a 
preferred direction at the xy orientation plane. The parameters and measure the 
substrate-fiuid interaction strength of each substrate. 

In our study, we will assume that the nematic director orientation is embedded in the xz 
plane. The formalism permits to define the orientation distribution of the molecules referred 
to a laboratory reference system, described by three order parameters, 

ri{z) = j d(psm9def{z,n)P2icose) (6) 
u{z) = j d4> sin 9def{z, n) sin 29 cos cf> (7) 
a{z) = J d(psin9d9f{z,n)sin^9cos2(p (8) 

The order also can be referred to an intrinsic reference system by three quantities: ip the 
tilt angle defined like the angle formed by the director with the z axis, U (uniaxial order 
parameter) the amount of order along the director and B (biaxial order parameter) which 
measures the amount of order along the perpendicular directions of the director. The nematic 
ordering is essentially uniaxial {B = 0) far from the substrate. The biaxiality induced by 
the surfaces was already discussed in a previous paper—. 
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These intrinsic order parameters can be obtained in terms of the set {77, cr, 1/} by the 
following expressions, 



tan ip{z) 



viz) - ^ + 



viz)-^ +iyizy 



1 2 




viz) - 



cr[z 



1 2 



+ uizY^ 



(9) 
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It is convenient to define the average density and order parameter for the confined system 
as: 



1 

H 



H 



piz)dz 



H 



U^jj I U{z)dz 



(12) 
(13) 



As we will also make reference to the semi-infinite problems, the adsorption in the order 
parameter is defined as: 



{U - Ub)dz 



(14) 



where f/f, is the value of the intrinsic nematic order parameter of the bulk phase. 

Numerical values for the potential parameters were taken as e^i = 1 (which sets the 
temperature scale), es/e^ — 0.847 and ec/^A — 0.75. The range parameters Aj are set, 
in units of d (throughout we choose this unit to set the length scale), to = 2,4,1.75, 
i — A,B,C respectively, and Xw = 1. In our calculations the temperature is fixed at the 
value T* = ksT/eA — 0.57, that corresponds to a typical situation in which vapor, isotropic 
and nematic phases can be observed. At this temperature, the isotropic-nematic coexistence 
occurs for a reduced chemical potential = /x/ZcsT = —3.918. 



III. THE SEMI-INFINITE CASE. 



In order to better understanding the behavior of the fluid inside the pore, we have first 
considered the semi-infinite cases of the fluid in presence of a single substrate that either 
favors homeotropic or planar alignment, i.e. the substrate-fluid interaction potentials are 
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given by Eq. (jU with e^y > or < 0, respectively. We are interested in the wetting 
behavior, so the system will be at bulk nematic-isotropic coexistence, unless stated otherwise. 
We have considered both situations in which the bulk fluid is either isotropic or nematic. 

First we consider the case of a nematic bulk phase in presence of the substrate. This 
situation was studied for substrates that favor homeotropic anchoring in Ref. We extend 
this analysis to substrates which favor planar alignment. For both negative and small positive 
values of e]y the fluid far from the substrate orientates parallel to the substrate. For the case 
in which < 0, biaxiality is developed in a microscopic layer close to the substrate, leading 
to the decay of the orientational order in that area. If the substrate favors a homeotropic 
alignment, we observe a microscopic isotropic-like layer close to the substrate, corresponding 
to a partial wetting situation of the substrate-nematic interface by the isotropic phase. As e]y 
is increased, the isotropic layer close to the substrate becomes macroscopic as the isotropic 
phase completely wets the nematic-substrate interface at e]y = . Our calculations showed 
that the wetting transition is continuous. The orientation far from the substrate remains 
planar in the complete wetting case. By further increasing e^, an anchoring transition occurs 
for ejy = ey[/, as the nematic orientation changes from planar to homeotropic far from the 
substrate, and the wetting becomes partial again since the isotropic-nematic interface favors 
a planar anchoring in our model. The anchoring transition is first order and continues in the 
single nematic phase region. However, the complete wetting with planar orientation state 
remains metastable up to e^. This fact will be important for the discussion of the confined 
cases. For the model parameters considered in this paper, e^} /^a = 0.470, ^w/^a = 0.527 
and e'^/eA = 0.72. 

We turn to the case in which there is an isotropic phase in bulk. For small values 
of \ely\, the inhomogeneities of the density and order parameter profiles are restricted to 
a microscopic layer close to the substrate. Although the fluid is isotropic far from the 
substrate, the anisotropic character of the substrate-fluid potential induces a random planar 
ordering near the substrate, i.e. the particles orientate parallel to the substrate but without 
any preferred direction in that orientation plane [i] < 0, a = u = 0). Such a condition can be 
seen as an extreme biaxial case, where B = 2U, i.e. the highest eigenvalue of the orientational 
order parameter tensor is twofold degenerate. The layer width remains microscopically finite 
as the NI coexistence is approached, corresponding to a partial wetting situation. For large 
enough values of |e^|, we observe complete wetting situations. For negative values of and 
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at the NI coexistence, there is a sudden change at ejy = e^y from the previous surface state 
to a complete wetting state of the isotropic-substrate interface by the nematic phase. This 
corresponds to a first order wetting transition. For our model parameters, j^A = —0.359. 
Due to its first order character, there is a prewetting first-order transition for ew < (see 
Fig. P). At this transition the particles at the first layer orientates along some direction 
parallel to the substrate, breaking the rotational symmetry of the state (see insets of Fig. 
2]for the order parameter profiles). Consequently, the prewetting transition can be seen as 
a two-dimensional isotropic to nematic transition at the fiuid layer close to the substrate. 
We will denote I and to the surface states with the random planar and nematic layer 
close to the wall, respectively. For the state, there is also biaxiality near the substrate, 
as occurred when there was a nematic fiuid in bulk. A similar behavior was found for the 
Zwanzig model and hard sphero cylinder fiuidslS, although their "prewetting" transition is 
continuous (second order for the Zwanzig model and it is expected to be Kosterlitz-Thouless- 
like for the sphero cylinder system). The prewetting surface nematization transition shifts 
towards lower chemical potentials as ew is increased, and eventually disappears at a critical 
point. 

For large positive values of e^, i.e. e^/e^ > e^/e^ = 0.691, we also found that the 
nematic phase completely wets the isotropic-substrate interface. We have checked indirectly 
this by studying the numerically obtained profiles that minimize our free energy functional. 
The order parameter profiles that we obtain in our numerical minimization are inhomoge- 
neous along the wetting layer (see Fig. Ej). This fact is due to the mismatch of anchoring 
conditions at the interfaces. We first note that the intrinsic order parameters are homoge- 
neous along the wetting layer except close to the interfaces (see inset in Fig. E}. So the origin 
of the inhomogeneity of the extrinsic order parameters is the change of the orientation of 
the nematic director field along the wetting layer. As previously mentioned, the molecules 
align parallel to the interface between a nematic and isotropic phase. On the other hand, the 
molecules anchor homeotropically close to the substrate due to the strong anchoring cou- 
pling between the substrate and the fiuid. As the nematic layer width diverges, it is possible 
to fulfill both anchoring conditions by slowly varying the nematic director from homeotropic 
to planar anchoring. However, any numerical calculation requires minimization in a finite 
range R for the z coordinate, imposing bulk isotropic conditions for z > R. In this situation, 
the nematic director profile that corresponds to the minimum of the grand-canonical free 
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energy is expected to be linear, as it is indeed observed by numerical minimization (see 
inset of Fig. Hj). This fact induces a. 1/R dependence of the apparent equilibrium excess 
grand-canonical free energy (or surface tension E) due to the elastic deformation. Such de- 
cay is much slower than the usual exponential correction when elastic deformations are not 
involved, and we have to consider much larger values of R to obtain reliable values of the 
true equilibrium free energies. As R increases, the apparent surface tension decreases and it 
is always larger than the complete wetting state excess free energy, that it is the sum of the 
surface tension S^^v corresponding to the interface between substrate and a nematic phase 
oriented homeotropically far away, and the surface tension E^v/ between corresponding to 
the nematic-isotropic interface where the nematic phase is planar-anchored with respect to 
the interface. The latter surface tension is independent of the fluid-substrate interactions, 
and its value is (3T,ni(P 0.0259 for the NI coexistence at T* = 0.57. 

Fig. El summarizes the results of this Section. We represent the cosine of the contact 
angle 6' of a droplet of isotropic phase in coexistence with bulk the nematic phase on a 
planar substrate as a function of ew- The Young equation states that: 

r, '^wN ~ ^wl -s, 

cosO = (15) 

where E^^r, E^/ and E^v/ are the surface tensions corresponding to the substrate- nematic, 
substrate-isotropic and nematic-isotropic interfaces, respectively. For ew < (including 
negative values), the nematic director is parallel to the substrate far from the substrate. 
Otherwise, the nematic director is homeotropically anchored. The complete wetting region 
corresponds to cos^^ = 1. The curve meets the cos 6 = 1 line at ew = ^w "with no change 
in slope. This fact is a signature of the second-order character of the wetting transition. 
On the contrary, the curve departs from the cos^^ = 1 line at ew = ^w "^ith non-zero slope, 
indicating the first-order character of the dewetting transition. 

Fig. 121 also gives information about the wetting properties when there is an isotropic 
phase in bulk. By reversing the roles of the nematic and isotropic phases in Eq. ()15|) . it is 
straightforward to see that the contact angle 6' of a nematic droplet on a planar substrate 
is related to 6 via 6' = tc — 6. Consequently, when cos^ = —1 we observe complete wetting 
by nematic phase of the isotropic-substrate interface. This fact is clear for ew < < 0- 
However, for e > the apparent values of cos 6* are less than —1. As previously discussed, 
this is a finite-size effect of our numerical calculations and the predicted value is cos 6 = — 1. 



9 



Nevertheless, since the nematic layer remains finite as ew ^yv from below, the localization 
of the wetting transition is quite accurate. For both wetting transitions at and there 
is a change of slope of the cos 6' curve, that it is in agreement with their first-order character. 

IV. THE CONFINED CASE. 

We turn to the case of the fiuid confined in a hybrid slit pore, where each substrate favors 
an antagonist anchoring, i.e. homeotropic and planar, respectively. Under confinement, it 
is well known that the bulk transitions have lower- dimensional counterparts that are known 
as capillary transitions. The thermodynamic conditions (temperature, pressure, chemical 
potential) for such transitions are usually shifted with respect to the bulk coexistence values. 
Before considering our DFT results, we consider the predictions of the Kelvin equation, 
which gives the leading order for the shift on chemical potential A/i = yU^j — /x^j of the 
capillary nematic-isotropic transition /i^j- with respect to the bulk value /i^j at large pore 
widths H. Thermodynamic consider at ions^iSSjSi show that the chemical potential shift A/x 
of a capillary transition with respect to the bulk one between the a and {3 phases at a given 
temperature obeys asymptotically: 



In this expression, Tjap is the surface tension for an a — /3 interface, and pa and pp are the 
bulk densities of the a and /? phases, respectively. Finally, 6i and 6*2 are the contact angles 
corresponding to the droplets of a phase that nucleate on each substrate when the fiuid is 
in bulk in the (3 phase, respectively. In our case, we can identify a as the isotropic phase 
and /3 as the nematic phase. It is implicitely assumed in Eq. (fT^ that there is a partial 
wetting situation for both substrates in semi-infinite geometries. However, if either there is 
complete wetting only at one substrate or in both but by the same phase {6i = 6*2 = or 
tt), the effective pore width in Eq. (fTBj) is reduced to H — nl, where k is a constant that 
depends on the range of the wall-fluid interactions and I is the sum of the widths of the 
adsorbed wetting layers. Nevertheless, this modiflcation does not change qualitatively the 
results. Its only effect is to observe the behavior predicted by the Kelvin equation at larger 
values of the pore width H than for a partial wetting situation. 

A different situation arises when both substrates are completely wet by different phases at 




(16) 
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coexistence, i.e. cos cos ^2 = —1- Studies for symmetrically opposing surface fields Ising 
systems^"^ and Landau-de Gennes models of liquid crystals^! show that in this situation 
there is no capillary transition. Instead, an interface-like state appears, where half the pore 
is filled with an a-like layer, and the rest by /5-like phase. The position of the interface 
between both layers is controlled by the deviation of chemical potential with respect to 
the coexistence value, being possible to go from an a-like capillary pure phase to a (3-like 
one without undergoing any phase transition. The phase transition is only recovered in 
the bulk limit H = 00. When the opposing surface fields are not symmetrical (at least for 
near-symmetrical), it is expected to observe a similar behavioi^. 

This analysis allows us to envisage different scenarios for the confined case and large H. 
For simplicity, we will only consider the case = = ew in Eqs. (jH) and (0). We 
must note that this choice does not correspond to the symmetrical opposing fields, since 
the effect of each substrate is qualitatively different. For large H and ew/^A < 0.48 and 
^w/^A > 0.53 we expect to have a nematic-isotropic capillary transition. The transition point 
will shifted towards lower or higher chemical potentials than the corresponding to the bulk 
transition if cos0i+cos^2 is negative or positive, respectively. This state-dependent capillary 
transition shift, which does not occur for simple fluids, has been observed in our model under 
confinement between symmetrical substrates that favor homeotropic anchoringi^, and also 
in confined spherocylinder model with the Onsager approximation for a symmetrical slit 
pore^S. On the other hand, if 0.48 < ew/^A < 0.53, no capillary transition is expected to be 
obtained for any pore width. 

We present now our DFT calculations. For the sake of clarity, we separate the results 
obtained for ew < > ^w- 

A. The case ew < 

We performed our first DFT calculations for ew/^A = 0.30, which correspond to partial 
wetting situation in the semi-infinite geometry for both substrates. Fig. |3]shows the behavior 
of the averaged intrinsic order parameter U as the chemical potential is varied. We see that 
there is a transition from an isotropic-like state (lower branch) to a nematic-like state (upper 
branch). The order of the transition is first-order, and it is found as the coincidence point 
of the grand-canonical free energies of the nematic-like and isotropic-like states. We have 
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shown also the metastable continuations of each branch: nematic-hke states on the left of the 
capillary NI transition and isotropic-like states on the right of the transition. Examination 
of the order parameter profiles show that the director field in the nematic-like state is 
oriented parallel throughout the pore to the substrate with an intrinsic order parameter 
similar to the bulk value, except in thin layers around the substrates, where the density and 
order parameter profiles are similar to the ones corresponding to the semi-infinite cases. In 
particular a very thin layer close to the z = substrate is homeotropically anchored. This 
state corresponds to the S state in Ref. |2l| and we will denote it by A^'^. On the other hand, 
the isotropic-like state (which we will denote as /) develops microscopic homeotropically 
and random-planar layers close to the z = and z = H substrates, respectively, and 
the orientational ordering decays to zero in the middle region of the pore. We have also 
compared our DFT calculations with the Kelvin law prediction Eq. (fT^ (see inset in Fig. 
H}. Although the transition value departs from this prediction for small pore widths, for 
H/ d > 20 the agreement between our numerical estimations and the Kelvin predicted value is 
quite reasonable. Note that the capillary NI transition is shifted to larger chemical potentials 
than the bulk NI transition value. However, as ew increases, this trend is reversed (the Kelvin 
law predicts the turning point to occur at e^/^A = 0.329). 

As ew exceeds e|^, wetting in the z = H wall begins to affect the pore adsorption. For 
small values of H a. single transition between an isotropic-like / state to a nematic-like 
state occurs as for ew < ^w- However, as H increases the order of the "nematic-hke" 
state decreases and the chemical potential for which this transition occurs converges to 
the I — transition value for the z = H wall. The order-parameter profiles make clear 
that the "nematic-like" state corresponds to the formation of a nematic layer close to the 
z = H substrate followed by an isotropic ordering along the pore. By analogy with the 
semi-infinite situation, we call to the confined state after this transition, which we can 
see as a reminiscence of the I — transition in the semi-infinite phenomenology. 

As H increases, in addition to the I — transition, a new phase transition between 
the to a proper nematic-like state A^'^ (i.e. planar- anchored along the pore) is observed. 
This transition can be regarded as the true capillary NI transition. Fig. El shows our DFT 
calculations for H = 20d and different values of ew It becomes clear that the difference 
between the and the A^'^ states at the transition decreases as ew increases, and eventually 
the NI capillary transition disappears in a critical point (see inset of Fig. Ej). The value of 
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ew at which the NI capillary transition disappears is lower than , i.e. the confinement 
distahilizes the two-phase region with respect to the one-phase region. We have also studied 
the effect of the pore width in the NI capillary transition. We set ew/^A = 0.465, slightly 
below . Fig. shows that the capillary NI transition does not exist for the smallest 
values of H and reappears for H > 50d. 

The disappearance of the NI capillar transition shown by our model is not related to the 
critical points observed for small values of H in other model o^^i^" . Actually, those cases 
correspond to symmetrical slit-pores and consequently their phenomenology is not related 
to wetting properties of the substrates. In our model the NI capillary transition occurs at 
least for H > 20d when the substrates at 2; = and z = H are the samoi^. However we 
cannot rule out the possibility of critical points for very small H, but we will not consider 
them since the predictions of our model are not reliable in that regime. 

Our DFT results about the existence of the NI capillary transition show the opposite 
trend to the one predicted by a Landau-de Gennes model in Ref. 13, which predicted 
that the confinement stabilized the two-phase region with respect to the one-phase region. 
Moreover for large values of H the boundary for the non-existence region of the capillary 
transition was a triple point rather than a critical point. Actually the scenario shown by 
our DFT calculations is consistent with the phenomenology that presents a confined Ising 
model with opposing surface fields when the wetting transition is critical^. This analogy 
can be explained heuristically. For large H and ew close to e^, the typical configuration 
close to bulk coexistence corresponds to an interfacial state formed by an isotropic layer 
around the z = and a nematic layer around the z = H wall. In first approximation the NI 
interface position is given by the minimum binding potential which is the superposition of the 
contributions from the two semi-infinite systems^. As there is a complete wetting situation 
at the z = H substrate, we can consider its binding potential as a exponentially decreasing 
function of the interfacial position with respect to the z = H wall cexp[—{H — 
where c > and is the nematic phase correlation length at bulk NI coexistence. Note 
that the local minimum corresponding to the bounded state, if any, plays no role in our 
argument as soon as we are far enough from the first-order wetting transition corresponding 
to the semi-infinite case with the z = H wall. On the other hand, the binding potential 
corresponding to the z = substrate takes the usual form — aexp(— //^t) -|- bexp{—2l/C^h), 
where a oc — ew > 0, b > and is the isotropic phase correlation length at bulk 
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NI coexistence. So, the complete binding potential for the NI interfase, up to irrelevant 
constants, is: 

W{1) = h{H -I)- ae'~b + he'~b + ce ^ (17) 
where h oc AyU. A critical point must verify: 

a _J_ 2h -iL c 
h = -e - —e + -e (18) 

kb kb kb 

a L 46 _2i c -4?^ , , 

= + + (19) 

Instead of solving simultaneously Eqs. ()18j) and we will study the spinodal line, which 
is the solution of Eq. We consider solutions of the form I = xH + Al, where < x < 1 

and Al ~ ^b,^b- Substituting this ansatz in Eq. (fT^ . we obtain 

= - ae «i> e «b + 46e e 



/ f2 (1-X)H \ Ai 

In order to have two solutions for the Eq. ()20|). we need all the factors in brackets in 
Eq. to be of the same order. Consequently we obtain that x = C^b/i'^Cb + ^b) and 

a = a' exp[—H/{2C,l + ^fe)], where a' b c. The value of a' for which the critical point 
occurs corresponds to the case where the rescaled equation: 

Ai 2Ai £2 ^ 

= -a'e ib+Abe + c-^e^b (21) 

[Q 

has an unique real solution. For simplicity, we assume that C,b = ^b^ but similar results are 
found when both correlation lengths are different. We refrain to give explicit solutions and 
we just mention that indeed a'/c ~ 1 and the value of Al at the critical point is Ale ~ ^& 
for reasonable values of b and c. Consequently, the interfacial position at the critical point 
Ic ~ H/3 and the NI capillary transition shift — eiy ~ exp(— ///S^t) and from Eq. ()18p 
Afi ~ exp{-2H/3^b)- 

A numerical comparison of these predictions with our DFT calculations, although it would 
be very interesting, is extremely difficult due to the limitations of our numerical procedure 
and it will not be carried out in this work. However, at least qualitatively these predictions 
have been confirmed by the DFT calculations even for not so large values of H. 
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B. The case ew > 

Now we turn to study the confined system for ew > ^w- ^^^^ considered tlie case 
^w/^A = 0.70 and H = 20d. Different confined states are obtained as tfie cliemical potential 
increases (see Fig. [7j). The first cases Figs. [7|^a) and[7|^b) clearly correspond to a vapor-like 
and isotropic-like states, and we will denote them as V and J, respectively. There is in 
both cases a homeotropically-anchored nematic layer close to the substrate at z = and a 
random planar layer close to the substrate at z = H. As the chemical potential increases, 
the random planar layer ai z = H undergoes a transition to a planar biaxial state. Again 
this first-order transition is a reminiscence of the I — surface transition corresponding to 
the substrate at z = H, and we will denote both confined states / and by analogy with 
the semi-infinite phenomenology. As the chemical potential is further increased, the nematic 
layer increases and a state in which an isotropic and a nematic layer coexist is obtained (Fig. 
|71^c)). There is no phase transition between the isotropic-like and the nematic-like 
states, but a smooth transformation similar to the expected one when both substrates are 
completely wet by different bulk phases at coexistence. Finally, there is a transition to a 



new nematic-like state (the L state in Ref. |2l|, that we will denote as A^^), characterized 



by an almost linear distorsion from homeotropic to planar configuration on the tilt angle 
through the pore, although the intrinsic order parameter takes essentially the bulk value 
except close to the substrates (see Fig. El^d)). This transition between the A^^ and A^'^ state 
is first order and converges to the anchoring transition that is observed in the semi-infinite 
case with a substrate that favors a homeotropic anchoring as H ^ oc^^. A more detailed 
discussion about this transition can be found in Ref. 0. 

For larger values of H, the capillary NI transition appears as a transition between an 
state and a A^^ state for chemical potentials smaller than the bulk NI transition value. As 
ew approaches e^, the range of values of H for which capillary transition does not exist is 
wider. For ew/^A = 0.53 (close to the dewetting transition point) the capillary NI transition 
does not appear until H/d > 150. Fig. |H1 shows the continuous change of the nematic layer 
as the chemical potential increases for H = 20d, which is analogous to the behavior observed 
for ew/^A = 0.70. Fig. El shows how the curves of the averaged intrinsic order parameter as 
a function of the chemical potential change for increasing pore widths. For the narrowest 
pore {H = lOd), a single transition is observed between an isotropic state to an state with 
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some nematic ordering. However, as H is increased, it becomes clear that this transition is 
again a reminiscence of the I — transition in the semi-infinite case. The locahzation of 
this transition depends weakly on the pore width H. For H > 13d, the A^^ — A^"^ transition 
appears at chemical potentials larger than the bulk NI coexistence value. This transition 
does not appear for H = lOd, which is a clear indication of the existence of a critical point 
for the A^^ — A^"^ transition at small values of H. As H increases, the chemical potential 
value for A^^ — A^"^ transition shifts to lower values. Finally, we observe that the tJ curve 
becomes increasingly steep at the bulk NI transition as — oo. This fact is a consequence 
of the formation of interfacial-like states (as Fig. Efc)) in the proximities of the bulk NI 
coexistence chemical potential. 

The non-existence of a capillary NI transition in the situation described above is unex- 
pected. First, the finite-width effects are much more pronounced for ew > than below 
e^. As a consequence, we observe significant departures from the bulk behavior for wider 
pores and in a larger range of ew than when we approach the complete wetting situation 
from below (i.e. ew < ^w)- This fact suggests the emergence of a new length scale in the 
system. On the other hand, as the dewetting transition is first order, we could naively expect 
that under confinement and large enough H the capillary transition disappearance point is 
shifted to lower values of ew^^ (recall that the role of ew or the temperature is reversed in 
the dewetting transition with respect to the wetting transition). To understand this situation 
we recall that the confinement in a hybrid cell distabilizes the A^^ with respect to the A^'^. 
A simple phenomenological argument can be used to explain this fact^^. For large enough 
values of H, the excess free-energy per surface area of a A^'^ state is given by the sum of the 
surface free energies between each substrate and a nematic fiuid which is planar-anchored 
with respect to them. The finite-width corrections are expected to decay exponentially. On 
the other hand, the large H behavior free energy of the A^^ state is the sum of the sur- 
face tensions of the interface between the z = substrate and a homeotropically-anchored 
nematic fiuid and the corresponding between the planar-anchored nematic fiuid and the 
z = H substrate. For this state the finite-width correction is dominated by the elastic de- 
formations contribution and decays as K^tt'^/SH a.s H oo, where is a Frank-Oseen 
elastic constant^. So, for a given temperature, the A^'^ — A^^ transition is always shifted 
towards higher chemical potentials than the anchoring transition, and the convergence to 
it for large values of H is quite slow. Although this argument implicitely assumes that the 



16 



inhomogeneities are restricted to microscopic regions around each substrate, it is also valid 
if we allow the A^"^ state to develop an isotropic layer close to the z = substrate (as occurs 
close to the bulk NI coexistence). Consequently, for e'^ < ew < moderate values of 

H, the most stable nematic-like state close to the bulk NI coexistence is , and therefore 
there is no capillary transition. As H increases, the state becomes more stable than the 
A^"^ state near the NI bulk transition, and eventually the capillary transition between a 
state and a A^-^ occurs. 

An interesting feature of the NI capillary transition in this regime is that it is linked to 
the purely surface transition — A^^. As explained before, for small values of H there is 
no NI capillary transition but if H is not too small there is a A^*^ — A^'^ transition. As H 
increases, the chemical potential corresponding to the A^'^ — A^'^ transition shifts to smaller 
values, and eventually crosses the value corresponding to the bulk NI transition (see Fig. 
ITUI) . When this happens, the transition is between a isotropic- like state to the nematic- 
like A^^ state: the capillary NI transition is recovered. Finally, the Kelvin relationship Eq. 
(fTH|l predicts that the capillary transition value of jSfi approaches the bulk NI value from 
below. The DFT results for large H support this prediction. So, at fixed temperature and 
chemical potential, and changing H, reentrance of the isotropic-like state is observed (see 

Fig. CUD. 

V. CONCLUSIONS 

In this paper we have reported a DFT study of the NI capillary transition under confine- 
ment in a slit-pore where each substrate favors a different anchoring condition. We found 
that there is an interplay between the phenomenology found in the semi-infinite problems 
with the phenomenology which appeared in the confined problem. We have focussed in 
the non-existence of the NI capillary transition as a result of both substrates to be wet by 
different phases. If we approach this region from below, we observe a distabilization of the 
two-phase region with respect to the one-phase region. Moreover, the capillary NI transition 
finishes in a critical point. These findings are different from the ones reported by Quintana 
et af'^ and this discrepancy can be traced to the different order of the wetting transition. 
When we approach the non-existence range for the NI capillary transition from above, a 
different mechanism is involved in the finite-size effects of the capillary NI transition. Al- 
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though the wetting transition which drives the disappearance of the capillary transition is 
first-order, the two-phase region is still distabilized with respect to the one-phase region. 
The explanation for this is the existence of two nematic-like locally stable phases with a 
relative stability dependent on the pore width. One of them corresponds to a complete 
wetting situation, so links with the isotropic- like states without undergoing a phase transi- 
tion. So, as H increases, the phase transition between two nematic-like states transforms 
continuously to the NI capillary transition. As a consequence, a reentrant behavior for the 
isotropic-like state is observed by varying the pore width and keeping the temperature and 
chemical potential constant. 

Although we have restricted our study to an isotherm case, we anticipate that we would 
find a similar behavior if the temperature is changed instead of the value of ey\r. On the 
other hand, the condition — is not so restrictive as soon as the wetting transition of 
the z — H wall precedes the complete wetting at the z — wall. If such a condition is not 
fulfilled, our conclusions may change. 
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FIG. 1: The adsorption T plotted vs. the chemical potential for a semi-infinite geometry where 
the fluid is isotropic in bulk and the substrate favors homeotropic anchoring. The dashed line 
corresponds to ew/^A = 0.3 and the continuous line to ew/^A = 0.66. The I — transition 
for the latter case is represented by the dotted line. The insets represent the reduced density 
p* = pd? and extrinsic order parameter proflles ry, a and v for different chemical potentials and 
^w/(-A = 0.66. 

FIG. 2: Plot of the reduced density pd^ and extrinsic order parameter profiles r/, a and for 
^w/^A = 0.70 and Pp = —3.918, obtained by numerical minimization for R = 435(i. Inset: the 
intrinsic order parameters U and corresponding to the same state. The intrinsic order parameter 
B is zero except close to z = R (not shown for clarity). 

FIG. 3: Plot of the cosine of the contact angle 9 of an isotropic phase droplet on the substrate- 
nematic interface as a function of the fluid-substrate interaction strength ew (thick continuous line) . 
The filled circles correspond to the numerical estimates for ew > (see text for explanation). 

FIG. 4: Plot of the averaged intrinsic order parameter throughout the pore U as a function of the 
reduced chemical potential Pp for ew/^A = 0.3 and H/d = 10 (continuous line) and H/d = 20 
(dashed line). The dotted lines show the location of the capillary NI transition. Inset: the diamonds 
represent our DFT values of the chemical potential for ew/^A = 0.3 at the capillary NI transition 
for H/d = 10, 20 and 35, and the continuous line is the Kelvin law prediction, Eq. (|16() . 

FIG. 5: Plot of the averaged intrinsic order parameter throughout the pore U as a function of 
the reduced chemical potential (3p for H/d = 20 and ew/^A = 0.35 (dashed line), ew/^A = 0.41 
(continuous line) and ew/^A = 0.45 (dot-dashed line). The dotted lines show the location of the 
I — transition (left) and capillary NI transitions (right). Inset: The symbols are the coexistence 
values of U at the / — and capillary NI transitions for different values of ew- The continuous 
lines are to guide the eye. 

FIG. 6: Plot of the averaged intrinsic order parameter throughout the pore U as a function of the 
reduced chemical potential j3p for ew /^A = 0.465 and H/d = 10, 20, 40, 50 and 60. 
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FIG. 7: Plot of the density (p* = pd^) and intrinsic uniaxial order parameter ([/) for the fluid 
confined in a pore of width H = 20d and ew/^A = 0.70 and chemical potentials: (a) P/i = —4.7, 
(b) = -4.45, (c) = -3.95 and = -3.7. 



FIG. 8: Density (p* = pd^) and intrinsic uniaxial order parameter profiles (U) for H = 20d and 
^w/^A = 0.53. The chemical potentials run from P/j, = —4.03 to P/j, = —3.64 with a step of 
A{pij,) = 0.01. 



FIG. 9: Plot of the averaged intrinsic order parameter U throughout the pore vs the reduced 
chemical potential Pp for ew/^A = 0.53 and H/d = 10, 13, 20, 30, 40 and 80. The discontinuities 
of U{(3p) correspond to the I — (left) and — (right) transitions. 



FIG. 10: Plot of the reduced chemical potential (5p corresponding to the — N^ transition (or the 
NI capillary transition, see text for explanation) as a function of the reduced inverse pore width 
d/H for values of ew/^A = 0.49 (circles), 0.53 (diamonds), 0.60 (triangles up) and 0.66 (triangles 
down). The continuous lines are only to guide the eye. 
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I. RODRIGUEZ-PONCE et al, FIGURE 1 
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I. RODRIGUEZ-PONCE et al, FIGURE 3 
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I. RODRIGUEZ-PONCE et al, FIGURE 4 
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I. RODRIGUEZ-PONCE et al, FIGURE 5 
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I. RODRIGUEZ-PONCE et al, FIGURE 6 
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I. RODRIGUEZ-PONCE et al, FIGURE 7 
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I. RODRIGUEZ-PONCE et al, FIGURE 8 
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I. RODRIGUEZ-PONCE et al, FIGURE 9 
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I. RODRIGUEZ-PONCE et al, FIGURE 10 
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